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Abstract 

Conformally invariant scalar waves in black hole spacetimes which are asymptotically anti-de Sit- 
ter are investigated. We consider both the (2 -|- l)-dimensional black hole and (3 -|- l)-dimensional 
Schwarzschild-anti-de Sitter spacetime as backgrounds. Analytical and numerical methods show that 
the waves decay exponentially in the (2 -|- 1) dimensional black hole background. However the falloff 
pattern of the conformal scalar waves in the Schwarzschild-anti-de Sitter background is generally neither 
exponential nor an inverse power rate, although the approximate falloff of the maximal peak is weakly 
exponential. We discuss the implications of these results for mass inflation. 
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1 Introduction 



It is well known that the maximally extended Reissner-Nordstrom spacetime can be imagined as a collection 
of different asymptotically flat universes connected by different charged black holes [0]. Except for the 
Schwarzschild solution, all the special solutions of the more general Kerr-Newman class of spacetimes can 
have two horizons, the inner and outer horizons. Nevertheless gravitational theorists find that these dual- 
horizon black holes are unphysical because causality can be violated inside the hole 0. Moreover any 
radiation (either electromagnetic or gravitational in nature) that goes into this kind of black hole will be 
indefinitely blue-shifted at the inner (or Cauchy) horizon This effect has caused some to expect that 
this null hypersurface acts like a barricade to the other universes in the maximally extended spacetime. 

This infinite blue-shift phenomenon at the Cauchy horizon was first discussed by Penrose in the late 
sixties 0. At that time, people believed that any small energy perturbation on these dual- horizon black 
holes would destroy the Cauchy horizon because the perturbation is indefinitely magnified there, causing an 
infinite spacetime curvature at the horizon. Thus the null hypersurface would become a spacelike curvature 
singularity and the gateway to other universes is sealed. Both numerical and analytic approaches |^ 
suggest that this null hypersurface is perturbatively unstable. However, although singularities are found at 
the Cauchy horizon, they are not spacetime curvature singularities at all. By taking the diverging stress- 
energy tensor into account, Hiscock showed that the perturbation only turns the horizon into a so-called 
whimper singularity: all curvature scalars are finite but a freely falling observer crossing the horizon measures 
an infinite energy density 0. This kind of singularity is too mild to seal off the passage to other universes, 
and so resolution of the issue would necessitate a study that did not rely on perturbation theory. 

Poisson and Israel made a breakthrough in this problem by showing that the Cauchy horizon can turn 
into a scalar spacetime curvature singularity Unlike Hiscock who considered a Reissner-Nordstrom black 
hole irradiated by a flux of incoming radiation, Poisson and Israel imposed both incoming and outgoing fluxes 
of radiation on a Reissner-Nordstrom background. The outgoing flux, even if negligibly small in quantity, 
makes the inner mass function of the black hole inflate without bound at the Cauchy horizon. More precisely, 
the inner mass of the black hole diverges at a rate of exp[Kv)/vP near the Cauchy horizon. The factor l/v^ 
comes from the decay rate of the scattered radiation tail 0, |^, 0, where p > and t> — > oo at the Cauchy 
horizon. Regardless of the values of p and the surface gravity k > 0, the mass parameter always grows 
although the exponential rate is attenuated by the decaying effect of the radiation tail. This phenomenon is 
called mass inflation and expected to seal the inner horizon because the diverging mass parameter induces a 
scalar curvature singularity at the horizon. This result implies that it is inappropriate to maximally extend 
any dual-horizon black holes beyond the Cauchy horizon because the spacetime is unstable against energy 
perturbations there. It is generally believed that such perturbations (in the form of gravitational radiation) 
always exist in more realistic black holes which does not have perfect spherical or axial symmetry. These are 
scattered around the black hole forming incoming and outgoing fluxes but they will eventually decay away 
as a tail of late time radiation ||^. In this way, mass inflation is expected to prevent violation of causality. 

In addition to the Reissner-Nordstrom black hole, the mass inflation phenomenon has been found to 
take place in other black hole configurations [0, |11], |T2|, |T3|, |T^, These configurations are in 1-1-1, 



2-1-1 and 3 -t- 1 dimensions as well as in asymptotically non-flat spacetimes. All of these calculations 
assumed the inverse power-law decay for late time radiation as an ansatz to obtain the inflating mass 
function near the Cauchy horizon. Since the mass parameter is attenuated by the decaying effect of the 
radiation tail, it is important to understand the behaviour of the radiative tail in a spacetime other than 
Reissner-Nordstrom class. Mellor and Moss have shown that the radiation from perturbations in a de Sitter 



background exponentially decreases jl^- Strictly speaking this result has nothing to do with late time 
falloff because the global geometry extends beyond the cosmological horizon. However it indicates that the 
radiative falloff behaviour is sensitive to the presence of the cosmological constant. More recent work by 



Ching et al. [113 demonstrated that under certain circumstances the tail can be something other than the 
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simple inverse power-law. 

In this paper we will study radiative falloff in spacetimes that are not asymptotically flat. We find that 
the inverse power-law is not universally true, and that in some asymptotically anti-de Sitter spacetimes, 
the late time tail decays exponentially. The asymptotically anti-de Sitter backgrounds we will study are the 
(2 -|- l)-dimensional black hole [T^ and Schwarzschild-anti-de Sitter spacetime. 

The outline of our paper is as follows. In section 2 we review the structure of the (D+l)-dimensional scalar 
wave equation in spherically symmetric spacetimes that are not necessarily asymptotically fiat, and discuss 
our numerical approach towards solving it. In section 3 we verify that our numerical approach correctly 
reproduces the power-law falloff in asymptotically flat spacetimes, and we cross-check this analytically. In 
the next two sections we study the falloff behaviour in the (2 + l)-dimensional (or 3D) black hole background 
and in (3 + l)-dimensional Schwarzschild-anti-de Sitter spacetime. Concluding remarks and an appendix 
round out our work. 



2 The Wave Equation in D -\- 1 Dimensions 

We shall study scalar waves in different dimensions, since wave equations for higher-spin flelds are of a 
qualitatively similar structure [^. The (conformally coupled) scalar wave equation in D + 1 dimensions is 

V^^ = ^R^, (1) 

where ^ is an arbitrary constant. If ^ = (Z^ — 1)/{4D), this equation is conformally invariant. We simplify 
the problem by considering only static, spherically symmetric {D + l)-dimensional spacetimes with metric 

ds' = -N{r)df + ^+r'dnn-i\ (2) 

JM[r) 

where N{r) is the lapse function and dVlD^i^ is the metric of a (D — l)-dimensional unit sphere. Generalizing 
Zerilli's separability technique |20[ to D dimensions, we assume 



^ = r(i-^)/2^(t,r)F,^ . (3) 
The functions are the D-dimensional spherical harmonics which satisfy the equation 



D 



-lil + D-2)Y,^ . (4) 



The product -l{l + D — 2) is the eigenvalue of the operator which is the angular derivative operator. It 
is straightforward to show that equation ([^) gives 

-duip{t,r) + N{r)dr[N{r)drip{t,r)]- N{r)Ve{r)ilj{t,r) = 0. (5) 

The function Ve{r) is deflned as 

,,,,.) . ,R,£^ ^^(r) ^ - 'I 'f - ^' mr) _ 1 1 , ^ ^ - 3) ^ - 1) . (6, 

2r dr 4r^ 4r^ 

One can rewrite the wave equation (^) as 

-dtMt,r) + N{r)C[^{t,r)] = 0, (7) 
with the help of a spatial differential operator 

C = dr[N{r)dr]-V,{r) . (8) 
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Alternatively, if we introduce 

dr 

X 



N{r) 



(9) 



then d^) can be written as 

dtt^l^{t,r{x))-d:,:,ilj{t,r{x)) + V{r{x))ilj{t,r{x)) = 0. (10) 



or as 



du^ijj{u,v) = -^N{r{u,v))Ve{r{u,v))ijj{u,v) (11) 

using null co-ordinates u = t — x and v = t + x. The function V{r) (defined as V{r) = N{r) K(r)) plays 
the role of a potential barrier which is induced from the background spacetime geometry. Although the 
potential when written in terms of the tortoise co-ordinate x can be very complicated, equation (|TUp has the 
familiar form of a potential scattering problem. 

Equation (|l^) can be integrated numerically in a straightforward fashion by using finite difference method. 
First of all the D'Alembert operator du — d^x can be discretized as 

ipit- At,x) -2ilj{t,x) +ilj{t + At,x) ilj{t,x - Ax) -2ilj{t,x) +ilj{t,x + Ax) ^, . 2\ ^/a 2^ 



At" Ax' 



0{At') + 0{Ax') (12) 



using Taylor's theorem. In order to formulate a well-posed Cauchy problem we need to include the initial 
conditions, which for simplicity we choose to be 

^(t = 0,a;) = and dt^{t = {),x) = u{x) . (13) 

Because the field ijj is initially zero, its subsequent evolution is solely the result of the initial impulse of the 
field dtip. Discretizing the second condition in ( |l3l) yields 



= .(x) + 0(At^), (14) 

where we employ a Gaussian distribution with finite support for u{x). We further define 

ip{mAt,nAx) = i)m,n , (15) 
V{nAx) = Vn, (16) 
u{n Ax) = Un , (17) 

where the mesh size has to satisfy the condition Ax > At so that the numerical rate of propagation of data 
is greater than its analytical counterpart. The discretization of the Cauchy problem above then implies 

1p-l^n = -AtUn, (18) 
i^O,n = 0, (19) 
At' 



m+l,n 



2-2 —7 - At' Vn 
Ax' 



At' 

'4^m,n ~ V^m-l.n + ~: 2 [4'm,n-l + 4'm,n+l ] ■ (20) 

Ax 



As a result, we can follow the evolution of the field ip starting from the initial data given at time t = 0. 

In the case where the black hole geometry is asymptotically fiat, the tortoise coordinate x goes from 
negative infinity to positive infinity. Therefore our Cauchy problem is similar to the infinite string problem 
in which the initial data propagates towards left and right indefinitely. The initial data no longer enjoys this 
privilege when the background is asymptotically anti-de Sitter because the tortoise coordinate goes from 
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minus infinity to zero only. In other words, the right-propagating data cannot travel in this direction forever. 
Analogous to the semi-infinite vibrating string problem, boundary conditions at spatial infinity (i.e. x = 0) 
are needed in the asymptotically anti-de Sitter background in order to formulate the problem appropriately. 
Two types of boundary conditions that are widely used in anti-de Sitter backgrounds are the Dirichlet and 
Neumann conditions [^. In our case, the former reads 



^(t,x = 0) = and d^ilj{t,x = 0) = 1 (21) 

whilst the latter is simply 

tlj{t,x = 0) = 1 and d^tlj{t,x = 0) = 0. (22) 

We shall consider employing both of these boundary conditions at spatial infinity for our numerical compu- 
tations whenever the background geometry is asymptotically anti-de Sitter. 



3 Asymptotically Flat Backgrounds 

In this section we will review the behaviour of radiative falloff in asymptotically flat background spacetimes 
[^. We will present the results of the numerical calculation first. 

Figure |l] shows a sample inverse power-decay of a scalar wave in the Schwarzschild background spacetime 
of mass M. We solve the wave equation (p!0| ) numerically using the scheme discussed in the previous section. 
The compact initial Gaussian impulse is centered at a distance r = 10 M (or a; = 12.76 M) and for simplicity 
we choose the I = 1 spherical harmonic. Figure |I| shows how the magnitude of the scalar field at a distance 
r = 20 M (i.e. x = 24.40 M) evolves. Using linear regression, we find that the slope of the straight line on 
the graph is -5.026 in agreement with the analytic prediction of an inverse power-law falloff with exponent 
2/ + 3 i. 

Figures |^ to ^ show the same potential barrier V{x) that is responsible for this falloff behaviour. Figure 
1^ shows an exponential decrease of the left side of the potential function; this is a result of the fact that the 
event horizon is located at x = - oo. On the other hand, the power-law decrease on the right side of V{x) 
shown in figure ^ is a direct consequence of the use of an asymptotically flat background. If we compare 
figures || to § with figures |^ to |^ which represent the potential function in Schwarzschild-de Sitter (SdS) 
spacetime, we find that the right side of V{x) has different decaying behaviour even the overall appearance 
of V{x) on the linear graph is very similar. It is this difference that distinguishes the falloff behaviour in 
the two backgrounds 

Now let us consider a scalar wave in an asymptotically fiat {D + l)-dimensional background. For the 
remainder of this section we will restrict our attention to the case where the number of spatial dimensions D 
is odd. The motivation for this may be traced back to Huygen's Principle, which implies that in even spatial 
dimensions the scalar wave obeying the equation V^\E' = always develops a tail, regardless of whether or 
not the asymptotically fiat background is sourceless. Consequently identification of the tail part of the wave 
that is due to solely to backscattering becomes quite problematic when D is even. 

Inspired by the work of Ching et. al. [|1^, we consider a spherically symmetric static metric of the 
form (0) with lapse function 

f Inlrl )^ 

N{r) = — ■ (23) 

The constants a and f3 are integers, where /3 > but a > 0. The other constant m is a real number. When 
we have a = 1, /3 = and D = 3, this becomes a Schwarzschild background of mass m/2. Once we have 
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Schwarzschild data: M = 1 , 1 = 1 , Ro = 10M, RO = 20M 




Figure 1: The decay of a scalar wave in Schwarzschild background. Prior tot = 200 the decay is accompanied 
by 'ringing' of the quasi-normal modes, after which the falloff rate is that of an inverse power-law. 



Schwarzschild data: M = 1 , 1 = 1 




Figure 2: Potential barrier for the Schwarzschild background 
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Schwarzschild data: M = 1 , 1 = 1 
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Figure 3: Exponential decrease on the left side of the barrier in V{x) 



Schwarzschild data: M = 1 , 1 = 1 




Figure 4: Inverse power decrease on the right side of the barrier of V{x) 
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SdS data: Lambda = 0.0001 , Rb = 2 (M = 0.999867), I = 1 , Rc = 86.09821 Rb 




X 



Figure 5: Potential barrier for the Schwarzschild-de Sitter background 



SdS data: Lambda = 0.0001 , Rb = 2 (M = 0.999867), I = 1 , Rc = 86.09821 Rb 
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Figure 6: Exponential decreasing nature on the left side of V{x) in SdS background 
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SdS data: Lambda = 0.0001 , Rb = 2 (M = 0.999867), I = 1 , Rc = 86.09821 Rb 
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Figure 7: Exponential decrease of the right side of V{x) in SdS background 



the lapse function, we can compute the potential function 

l-a){D 



VJr) 



(2l + D-3)(2l + D-l) AHD 
h m 



a] 



{D-l){D-3-2a) (ln|r| 



+ m P 



4^2 

2^{2D-3-2a)-D + l (ln|r| 



N/3-1 



M+2 



+ m/5(/5- l)e 



Inlrl 



(24) 



from 



We first find the static solution ips{r) of wave equation (|), by obtaining the solution of the equation 
C[i/js{r) ] = 0. It is straightforward to show that ipsi^) has the following form: 



ips{r) 



j=0 



r3 < 



7+1 ^ "^j (^) 



r3 < 



(25) 



where 7 = / + (-D — 3)/2. Notice that 7 is an integer when the spatial dimension D is odd. Except when 
D = 1 this integer is always positive. Hence it is the first sum of the solution that is physically relevant 
since it vanishes for large r, and we choose this as the static solution. For the remainder of this section, 
we assume that D > 3 so that this choice of ips{f) is valid. The coefficients ao('") and co(r) are arbitrary 
constants but the other coefficients are all polynomial in ln|r|. The generating equations for aj{r) and Cj{r) 
are given in Appendix 0. 

We follow the approach in refs. P, E31 and let 



ijjj 



i=0 



l^-^\u) + {-lYf-^\v) 



(26) 



be the form of the initial wave, i.e. the wave emitted by a star at the onset of gravitational collapse. In other 
words, this is the time when t <^ r. The functions g{u) and f{v) are as yet unknown. The term g^~'^\u) 
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represents i integrations of the function g{u) with respect to u; similarly for Using 
([IT]) becomes 

2 ar -I 
This equation has a set of solutions 



equation 



Bn(r) 



\N{r) ^B,{r)-\ 
2 ar 2 



K(r) B,{r) dr 



0,1,2, 



(27) 

(28) 
(29) 



where we have set Bo{r) = constant = 1 without loss of generality. The pair of equations above allow us to 
generate Bi{r) hierarchically in a straightforward manner. 

We can split each Bi{r) into two parts, denoted by Bf{r) and Bj{r). Bf{r) is defined as the m- 
independent portion of Bi{r), while Bj{r) is the rest which is m-dependent. Physically, the part Bf{r) 
represents the wave on the lightcone because it is the part that would be generated if the background were 
fiat (m = 0). Price referred to this part as the primary wave which depends only on the mode of the spherical 
harmonics. The other part Bj{r) is called the tail of the wave because it is created by the presence of the 



spacetime curvature and is off the lightcone due to scattering. Given equations 
that the primary part of -Bj(r) is simply 



and (By), one can show 



Bfir) 



r(7 + 1 + i) 

2H\T{-i + 1- i)r^ 



(30) 



When D is odd, 7 is an integer and the sequence { Bf (r) }j=o truncates. However if D is even, the sequence 
does not terminate and the primary part of the initial wave ipi has infinitely many terms. Indeed when this 
is the case, it is inappropriate to call Bf{r) the primary part because there are tails present. The approach 
breaks down because the primary and tail parts of the wave become indistinguishable. This is the reason 
we restrict ourselves to odd D in this section as mentioned earlier. 

Let us now consider the tail part of Bi{r). Unlike the primary part, Bj{r) has no simple solution. 
Fortunately one can always generate Bi{r) recursively. Note that Bf{r) is of order 0{ (ln|r| Y /r*+° ). 

Since 7 is an integer in odd spatial dimensions, it is convenient to define the functions G{u) and F{v) as 



G{u) = g^-^\u) 



and 



(31) 



such that the initial wave can be expressed as ipi = J2 Biir) Tp \u, v) where 



■y—i 



(32) 



When the background is fiat (m = 0), only the first 7 + 1 terms of Bf{r) survive. We expect that any 
outgoing radiation will propagate to spatial infinity without any scattering in a flat background. Therefore 
we must have 



_p(7)( 







(33) 



because F{v) represents the scattered infalling radiation and the Bf (r) are non-zero only for < i < 7. In 
other words, we insist that F{v) = 0. 
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For concreteness, we suppose the scalar wave starts leaving the star at retarded time u = Uq; that is to 
say the scalar field is zero before u = Uq. Continuity then implies 

GW(f/o) = G^-'-^XUo) = ■■■ = G'-^\Uo) = G^^\Uo) = (34) 

G^-'\u) = / G^'\OdC, G^-'\u) = / G(-i)(C)rfC, ••• (35) 
juq juq 



We also assume that there is essentially no emission from the star after u = Ui > Uq (due to gravitational 
red shift). When the star possesses no initial static moment at the onset of collapse, this cutoff condition is 
translated into 

GW(m) = G^^-^\u) = ■■■ = gW(w) = G^'^Ku) = 0, yu>Ui>Uo (36) 

because only the terms G{u) to G^'^\u) which represent the primary wave correspond to the emission from 
the star. This condition implies 

G^-^\u) = rG^^\C)dC = r'G^°\C)dC = constant, W u > . (37) 
Therefore if there is no static moment before u = Uq, the scalar wave after time u = Ui is 

oo 

'^u>u^ = E B,ir)Ti^~\u,v) , (38) 

4=7+1 

where Tf^{u, v) is only a constant and Bi(r) = Bf{r) is of order 0{ (ln|r| )^ y go the first term of 
this solution which is time-independent is of order 0( (ln|r| Y /r'''+"+^ ). Physically, this says the tail of the 
perturbation persists after the time u = Ui. The primary part of the wave has been "washed out" by the 
gravitational red shift before u = Ui. 

The situation is more complicated if the star carries a static moment at the onset of the collapse because 
the static wave ipsi^) can also be divided into primary and tail parts. In this case, we superimpose ipsif) 
of (^) and ipi{r) of (^) together to form a new initial scalar wave. We also suppose that only the tail of 
the perturbation persists after the retarded time u = Ui [^]. For the superimposed initial wave this cutoff 
condition requires 

G^''\u) = Vn>t/i>[/o (39) 

(27)! 

which yields 

G'W(m) = G^^-^\u) = ■■■ = G(^)(m) = 0. (40) 

The restriction on G{u) above ensures that the primary part of (^) can be cancelled by the m-independent 
part of ( P^D properly after the time Ui, leaving the combined ipi to be m-dependent only after Ui. As a 
result, the wave with initial static moment will become 

^u>v. = ^t"i^ + tBr{r)T,^-{u,v). (41) 

j=l «=7 

One can show that the first sum in the equation above is of order 0{ (ln|r| )^ /r"'+°' ) because the coefficient 
function ai(r) is a polynomial in ln|r| with degree p. In the second sum, the term Tj^{u, v) is a constant which 
equals to G{u), and B^{r) is of order 0{ (ln|r| Y /r'^"'"" ). Therefore instead of having a time independent 
part of order 0( ( ln|r| )'^ jcj^i+a+i ^^^q wave with initial static moment has an order 0( ( ln|r| )'^ /j"7+" ) only. 
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Now we turn our attention to the scalar wave at late time. By late time, we mean t ^ r. For the late 
time wave 7/'^, we introduce another ansatz 

00 

= J2C^{r)Tl{u,v) , Tl{u,v) ^ I^^\u) + {-ly H^^\v) . (42) 

i=0 



By substituting ip = ipi, equation ([TTD becomes 

1 1 f d ^ 

-Ar(r)£[C'o(r)]T°K.;) + -iV(r)5:|£[a+i(r)]-2-a(r)|71+iKt;) = 0. (43) 

This equation has another set of solutions 

C[C,{r)] = 0, (44) 

/:[a+i(r)] = 2-^a(r), 2 = 0,1,2,- •• (45) 
dr 

Unlike the case for the initial wave ipi, (for which the functions -Bj(r) can be calculated recursively using 
straightforward differentiation and integration), recursive generation of the functions Cj(r) involve inverting 
the differential operator C. The zero order equation has a solution 

^o(r) = r^+^E^- (46) 

The coefficient functions Cj(r) are those in equations (P^j). They can be calculated by using the generating 
equations in Appendix 0. The other inhomogeneous differential equation (|45|) has a solution of the form 



00 I 



Q(r) = r^-^^+^E^^ (47) 

where i = 0, 1, 2, ■ ■ ■ When i = 0, it is understood that c°(r) = Cj{r). The coefficients c^(r) are also given in 
Appendix ^ Since each coefficient CQ(r) is a constant instead of a polynomial in ln|r|, we can estimate the 
order of the late time wave tpi as 

00 

i^L = Y.O[r^''^^') Tl\u,v). (48) 

i=0 

Finally we match the late time solution to the initial solution at some transient period where u, v and 
r are of the same order. As the background is asymptotically flat, the tortoise coordinate x must have an 
order similar to that of r in this transient region. That is to say the orders of r and t are the same. The 
initial wave ipuyui ^^is period becomes 



Inlrl 



'4'u>Ui ~ O ( ^ yla+i ) ' without initial static moment 



Inlrl 



\/3 



ipu>Ui ~ ^ T — ) with initial static moment . 

In order to have consistent orders in the transient period, we must have 



TL'{t,t) ~ o(^ 



Inltl 



^/3 
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(49) 



if the star has no initial static moment; whilst the order of T^" must be 



T^{t,t) ~ o(^ 



Inltl 



(50) 



if the star has a static moment at the onset of collapse. In other words equations (|49| ) and ([50| ) give the 
falloff behaviour of the wave at late time. The inverse power falloff behaviour modified by a logarithmic 
term was first noted by Ching et al. E^. For a Schwarzschild background, we set D = 3, a = 1 and /5 = 



and have 7 = /, yielding the familiar power- law decay rate P, 23 



4 The 3D Black Hole Background 



The (2 + l)-dimensional black hole spacetime obtained by Banados Teitelboim and Zanelli |jT9[ is a spacetime 
which satisfies the vacuum Einstein equations with a negative cosmological constant (A < 0). Its metric is 



N{r) 



- N{r) dt^ + 



|A| r^-M + 



dr"^ 
iV(^ 

J2 



272 



dt + 



4^2 



(51) 
(52) 



and is actually anti-de Sitter spacetime with identifications. The constant M > is the quasi-local mass of 



the black hole and J is the angular momentum of the hole ||2^. For our purposes we shall regard the above 
metric as a spacetime which is asymptotically anti-de Sitter. 
When J = 0, the wave equation (0) gives (p!0D with 



3|Ak , M + 



4r" 



(53) 



In this case, the tortoise coordinate x is given by 

dr 1 



X = 



N{r) 



In 



AIM 



|A|r 



|A|r + VM 



(54) 



As r goes from yM/|A| to infinity, the tortoise coordinate has a range {-00 ^ 0). Thus we can write r in 
terms of 



l + exp(2 J|A|Mx 



(55) 



1^1 1 -exp(2^|A|Mx 

For conformal wave in 2 -|- 1 dimensions, the parameter ^ equals 1/8 and the potential barrier V becomes 

exp(2 A x) 



V{x) 



^1 + exp(2Ax) f 
where, for convenience, we have defined 

Vo = |A|(M + 4/2) > and 



Vo exp(2 A x) 



A 



Al M > . 



(56) 



(57) 



The procedure for finding a solution to the conformal scalar wave equation is as follows. We first write 
down the representation of the solution ip in terms of Green's functions. The rest of the problem then reduces 
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to that of looking for the correct Green's function. We will Fourier transform from the time domain to the 
frequency domain and obtain the Green's function in frequency space. Once we have the Green's function 
in frequency domain, an inverse Fourier transformation will yield the solution ip{t,x). The solution is exact 
for the approximate potential (|^). We will then use this exact solution and the Born approximation to 
compute the next correction of the solution for a more accurate description of the potential. A case when 
J 7^ will be studied as well. 

We first assume that there exists a Green's function 

G{x,C,t-r) = G{^,x-t-r) (58) 

which is zero when t < t. We define an operator D as 

D = du-d,, + V{x), (59) 

such that the Green's function with respect to this operator has the property 

DG{x,^;t-T) = [du-d,, + V{x)]G{x,^;t~T) = 6{t - r) 5{x - ■ (60) 

The inner product between Dip{t, x) and G{x, — t) gives 

ij{t,x) = f [G{x,^-t)dtmi) + mOdtG{x,i-t)]di 

J- oo 

+ / [G{x,i-t-r)d^,l,{T,0-i'{r,Od^G{x,^-t-T)f^^dT . (61) 



As a result, we have changed the question from looking for ip{t^x) to searching for an appropriate Green's 
function. 

We now carry out a Fourier transformation and define 

/oo 
G{x,i]t) eycp{iujt)dt . (62) 
oo 

Therefore equation (^) becomes 

DG{x,i-u) = [-u^-d^^ + V{x)]G{x,i-u) = 5{x-S) (63) 
if the Green's function G{x,C,',t) satisfies the conditions 

limG(x,^;t) = lim dtG{x,^;t) = 0. (64) 

t — ^oo t — >oo 

On physical grounds these assumptions are reasonable because any localized quantity is expected to be 
dispersed throughout the space by means of wave propagation. Mathematically, these assumptions are 
consistent with the Fourier transformability of the function G{x,C,',t) which must be absolutely integrable 
over 3fJ in order for G{x,C,', uj) to be well-defined. The Green's function in frequency space can be represented 
as 

G{x,^;uj) = (65) 

I W{L,;gJ) ^ ? , 

where the function W{uj]g, f) is the Wronskian of two linearly independent functions f{x]u) and g{x]Lj), 
that is 

W{uj;gJ) = g{x;uj)d^f{x;uj) - f{x;uj)d:,g{x;uj) . (66) 
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The functions f{x]uj) and g{x]uj) are two independent solutions of the equations 

D fix; 00) = Dg{x;uj) = 



(67) 



so that W{u!; g, f) is independent on x. 

Let ip{x;uj) be the Fourier transformation of the solution ip{t^x). We impose boundary conditions on 
the functions f{x;u) and g{x]uj) so that 



f{x;uj) oc ip{x;uj) and dxf{x;uj) oc dxip{x]uj) 

at the point x ^ - 00. At the other end (x = 0) we insist that 

g{x]u) oc iIj{x;u) and d^gix^u) oc dxip{x]Lj) . 



(68) 



(69) 

In equation (|6TD , after we have used the inverse Fourier transformation and interchanged the integrals with 
respect to du and dr, the representation of ipit^x) simply becomes 



ij{t,x) 



: G{x, 1) 9i^(0, + ^(0, dtG{x, ^; t) ] d^ 



(70) 



4.1 Solution for an Approximate Potential 



We will obtain the exact solution for the scalar wave corresponding to the approximate potential barrier 
( p6D in a static 3D black hole background. This entails finding the solution of the differential equation 
Dh{x;uj) = 0. This equation has two solutions [|l^, namely 



h{x; u) 
z{x) 



iJzix) 



Zq exp(Aa;) 



or K^{z{x)) 



A 



exp(Aa;) , 



(71) 
(72) 



where v = - iu / \, and Ii, and Ki, are the modified Bessel functions of complex order v . At the boundary 
X - 00 (the event horizon of the 3D black hole) we employ the condition 

f{x;u) = exp{-iujx) x^-00 (73) 

as usual |T^. Since this spacetime is asymptotically anti-de Sitter instead of flat, the boundary condition 



at spatial infinity is less trivial. Choosing Dirichlet conditions ^1 

g(x = 0; cu) = and 
the two functions f{x; u) and g{x; uj) read 



dxg{x;uj)\x=o 



f{x;uj) = Ufh{z{x)) , 

g{x-u) = Afg[I,{Zo)K,{z{x))~K,{Zo)lAz{x) 
where the normalization coefficients A// and Ag are given by 



^ff = tt) r(i + z/ 



and 



1 

A 



Now we put everything together and obtain 

/,( z{0 ) [ UZo) K,{ z{x) ) - K,{Z^) h{ z{x) 



G{x,i < x]uj) 



A 4(^0) 



(74) 



(75) 
(76) 



(77) 



(7J 
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The remainder of the problem is to bring it from the frequency domain back to the time domain. 
The inverse Fourier transformation is given by the equation 

1 roo ~ 

G{x,^;t) = — / G{x,^;iu) exp{-iLut)duj . (79) 
2 vr J-oo 

We evaluate this integral by analytically extending u to complex values and using Cauchy's residue theorem. 
The contour of integration is chosen to be a large, closed semi-circle on the lower half a;-plane, with center 
at = 0. By Jordan's lemma, the contribution from the arc of the circle goes to zero as the radius of the 
arc tends to infinity. The Green's function then becomes 

G{x,^;t) = i ^Res{G(x,^;cj) exp(-icjt) } , (80) 

where Res{ h(z) } denotes the residue of the function h{z) at a pole. In our case, the poles of the Green's 
function G{x,^] uj) come from the term Iu{Zo) only. This term causes trouble because it has infinitely many 
zeros in the lower half cu-plane. However z/ = -za'/A = Ois not a root of Iu{Zo) for positive definite Zq and 
so G{x,C,]t) decays to zero exponentially. The dominant exponential decaying rate is determined by the 
zero closest to the real axis on the lower half cj-plane. In general this root has a non-vanishing real part 
which contributes to the oscillatory nature of the falloff behaviour. As a result, the exponentially decaying 
Green's function G{x,^;t) implies that 'ilj{t,x) also decays in this manner. 

Had we chosen the Neumann condition instead of the Dirichlet condition for g{x,uj) |2l|], that is 

g{x = 0;uj) = 1 and d^g{x]uj)\^=o = 0, (81) 

the function g{x,uj) would simply become 

g{x-u) = Zo[IJ{Zo)K{z{x))-Kj{Zo)I,{z{x))] . (82) 

The change in Afg does not affect the result in any way because it is always cancelled in G when the Wronskian 
is taken into account. The change from Iu{Zo) to IJ{Zq) causes the singular term in G to become IJ{Zq) 
which also has infinitely many roots on the lower half tu-plane. Again, the origin is not a root. Therefore 
even if we employ the Neumann condition, the solution ip{t,x) also dies out to zero at an exponential rate. 

The next few graphs are the numerical results using the exponential potential. We used the same initial 
condition as in the numerical computation for the Schwarzschild case. Since the background is asymptotically 
anti-de Sitter, we integrate equation (p!0| ) numerically with the Dirichlet boundary condition at r = cxd (i.e. 
a; = 0.) Figures ^ and |^ illustrate the falloff behaviour of the wave using different exponential potential 
functions. In both graphs, the initial Gaussian impulse is located at twice the black hole radius {2Rb) 
and the observation is made at 47?^. The straight line asymptote of the ringing behaviour on the semilog 
graph corresponds to exponential falloff, numerically confirming that the wave exponentially decays in this 
asymptotically non-flat background. 

A direct check of our analytic prediction against the numerical computation is problematic because no 
analytic procedure exists for computing the (complex) values of u for which Iu{Zq) = for fixed Zq. Hence 
finding the root of /-j^/A(^o) which is closest to the real axis on the lower half cu-plane must be done 
numerically, by checking whether the root measured from the numerically generated graph (e.g. Figure ^ 
satisfies the equation /-^^^/^(^o) = for the relevant Zq (e.g. Zq = 3 for Figure |]). We have checked that 
the numerically determined roots satisfy this criterion for various values of Zq for all computations carried 
out in the (2 + l)-dimensional case. 

Figures |T^ and |TT| illustrate the numerical results using the approximate potential ( |5B| ) and the Neumann 
boundary condition at a; = 0. It is clear from the two graphs that the scalar wave also exponentially decays 
when the Neumann condition is used. 
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lambda = 0.02, Zo = 3, Ro = 2Rb, RO = 4Rb 

10" I , , , , , 




100 200 300 400 500 600 

t 



Figure 8: Exponential falloff for a conformal scalar wave using the approximate potential (^) and Diriclilet 
boundary condition 



lambda = 0.1 , Zo = 1 0, Ro = 2Rb, RO = 4Rb 
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Figure 9: Another semilog graph as in figure ^ but with different parameters 
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lambda = 0.02, Zo = 3, Ro = 2Rb, RO = 4Rb 
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Figure 10: Exponential falloff using the Neumann condition 
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Figure 11: The same as in figure ^ but with different parameters 
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4.2 Higher order correction 

The results of the previous subsection are based upon the use of the potential barrier V{x) = Vq exp(2 Ax) 
which is an approximation to 

V{x) = Vo ,,,\,2 - (83) 

( 1 + exp(2 Ax) ) 

Since A > but x < 0, exp(2 Ax) is always less than one. Therefore V{x) has a converging series expansion 

V{x) = ^0 [exp(2Ax) - 2 exp(4Ax) + 3 exp(6Ax) -4 exp(8Ax) + O(exp(10Ax))] . (84) 

As a result we may iteratively solve the wave equation, treating Vq exp(2Ax) as the lowest order term, 
-2 Vo exp(4 Ax) as the first correction to this approximation, and so on. 
For clarity, we define 

Voix) = roexp(2Ax), (85) 

Vi{x) = -2K) exp(4Ax) , (86) 

D = -uj^-d,, + Vo{x) + Vi{x), (87) 

Do = -uj^-d,, + Vo{x) . (88) 

Therefore the equation Df{x; u) = has a representation 

fixM = foix;u) + l^ W{u;;U.,U^) ^^iOm^dC (89) 

at X = - oo. The functions f/+(x) and U-{x) satisfy the equation DqU+{x) = DqU-{x) = 0. The other 
function /o(x; u) satisfies the same equation as ?7+(x) and U- (x) and also satisfies the boundary conditions 
at X ^ - oo. In other words, /o(x; u) is the solution (|75D in our case. Similarly the solution around x = is 

= ,o(x;.) - r "^^^llf ~V,T'''^^ VdO,iC,^)K, (90) 

Jx W[u!;U-,U+) 

where go{x; u) is just (^). We first compute the Wronskian W{uj] g, f) for /(x; u) and g{x; u) above. Since 
this Wronskian is x-independent, the simplest way to compute it is to evaluate the quantity at the point 
X = 0. It is not difficult to show that 

Wiu;gJ) = l^(c^;^7oJo) + W^(c^;^7oJ)(x)i:f^ = l^(cu;^?oJo)- W^(^;^7jo)(x)i:f^ . (91) 

Since g{x;uj) and f{x;uj) satisfy a differential equation which differs from that satisfied by go{x;uj) and 
/o(x; u), the Wronskians W{u!; go, /) and W{uj; g, fo) are functions of x in general. Therefore the correction 
Vi{x) in the potential barrier induces an extra (x-independent) term in the Wronskian W{uj; g, /). 
The first Born approximations for f{x;u}) and g{x\uo) read 

f{x-uj) ^ /i(x;^) = /o(x;cu) + ^ r^'^^[ir,(^(x)) /,(s)-/,(^(x)) ir,(s)] /o(s;u;)s=^c?s (92) 



Vo Jo 
and 

2 A^ f^o 

g{x]uj) ^ gi{x;u) = goix]Uj) - —— [K^{z{x)) I^{s) ~ I^{z{x)) K^{s)] go{s;uj) ds (93) 
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if we choose U- (x) = z{x) ) and U+{x) = z{x) ). The functions /o(s; u) and go{s] u) are defined via 
the equations 

fo{z{xy,uj) = fo{x;uj) = Ufh{z{x)) , (94) 

Uz{x)-uj) = g^{x-uj) = Ng[h{Z^)K,{z{x))-K,{Z^)h{z{x))] . (95) 

As a result, the Wronskian W{uj]g, f) can be approximated by 



W{u-gJ) ^ W{LO-go,h) + W{u:-g^,h){x)t^ = W{u-go,h)+ / go{x;uj) Ux-u)Vi{x) dx . (96) 







oo 



This Wronskian has roots on the lower half tu-plane. Thus there is quasi-normal ringing effect in the wave 
tail in this 3D black hole configuration with the approximate potential function V{x) ~ Vq{x) + V\{x). 
Moreover, since gi and /i have no cut, the falloff rate of the conformal wave is also exponential. 



The shape of the exact potential function V{x) is given on figure ^with parameters A = 0.02 and Zq = 3. 
Figures ^ and |n| are results of the numerical calculations using the exact 3D black hole potential (p6D. 
The methodology underlying these computations is identical to the previous case that uses the approximate 
exponential potential function. The parameters and boundary conditions used in this case are also the same 
as the cases in figures |^ and |^. Once again the numerical calculations show the exponential falloff of the 
scalar wave in the static 3D background. 




X 



Figure 12: The shape of the exact 3D black hole potential V{x) 



4.3 Spinning 3D Black Hole Background 

When J 7^ 0, i.e. the black hole rotates, we also find a late time exponential decay rate, as we shall now 
demonstrate. 



20 



BTZ data: lambda = 0.02, Zo = 3, Ro = 2Rb, RO = 4Rb 
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Figure 13: Exponential decay of the scalar wave ip in the static 3D black hole background 



BTZ data: lambda = 0.1, Zo = 10, Ro = 2Rb, RO = 4Rb 




Figure 14: Scalar wave decay as in fig 113] but with different parameters 
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If we assume that 



(97) 



that is to say there is no "spherical harmonic" component to the wave, the conformal scalar wave equation 
V^\E' = i?\l//8 will reduce to (|^). The tortoise coordinate x is defined as before but the potential barrier 
becomes 



Veixir)) 



3|A| 



— drNir) 
2 r 



N{r) 



In this spinning case, the black hole has two horizons R± which are given by the equation 

1 



R 



± 



2|A| 



M ± JM^ - |A| J2 



The lapse function N{r) can then be written in terms of -R+ and R- as 

Nir) = ^(r^-i?. 



As a result, the tortoise coordinate reads 

""^"^ ^ 2|A| (V-iL^ 
which has an inverse of the form 

R^ 



R^ In 



r-R^ 
r + R, 



R- In 



r-R. 
r + R. 



Y 



r{x) 
a 



1-^2 



2n 



n=0 



y 

A 



1 — exp(2 A x) 
1 + exp(2Ax) 
|A|i?+(l-a2) > 0. 



< 1 



The coefficients a„(cr) which read 

ao(o-) = 1 
25-17(72 + 3^4 



02 (cr) 



15(1 - a 



2^4 



ai((Tj 
03(0-) 



cr 



- a 



- a 



a 



2\2 



3(1 

1008 - 1039 or^ + 368 - 45 
315(1-^2)6 



are of order 0( (t2"' ). If a is small enough, we can employ the approximation 

1 _ Y 
rix) ^ R^ ■ 

If we substitute this approximation into the potential barrier V{x), we will obtain 
2 (i?+^-i?-^)(i?+^-4i?-2) exp(2Aa:) 



V{x) 



R, 



+ 12 A^ (2i?+2 - 3i?-' 



R- 



1 + exp(2Ax) f 
exp(4 Aa;) 



R+^ (1 + exp(2Ax) 



+ 80A^ 



-R-^ exp(6Ax) 



(98) 
(99) 

(100) 
(101) 

(102) 
(103) 

(104) 
(105) 

(106) 
(107) 

(108) 



R+^ (H-exp(2Ax)) 



6 • 



(109) 
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Figure 15: Potential function V{x) of two spinning 3D black hole backgrounds with (i) |A| J^/M^ = 10" ^ 
and (ii) |A| JVM2 = 0.8 

It is obvious that the first approximation of this potential has a form V{x) ~ Vq exp(2 Ax) while the second 
one is V{x) ~ Vq exp(2Ax) + Vi exp(4Ax). According to our previous analysis, we may conclude that the 
decay rate at late time is also exponential rather than a power-law. 

Figure |15| is the graph of the potential functions of this spinning case with zero angular harmonic. One 
can show that when the ratio |A| J^/M^ > 16/25, becomes negative for some r > R^. This behaviour is 
illustrated in figure where V{x) becomes negative when |x| is sufficiently large. In figure 0, we can see 
that the potential functions in both cases vanish at an exponential rate towards the event horizon. As usual, 
we put a Gaussian impulse at a distance twice of the outer horizon and the numerical response of the 
scalar waves over time at a distance 4 i?+ is shown in figure |1^. This figure shows an exponential decay of 
the wave although its appearance differs from the previous graphs because there are no angular harmonics 
(i.e. / = 0). Indeed if we set / = in the static 3D black hole case, we find that the falloff also looks like that 
in figure [l^ because the pole that is closest to the real cu-axis for the Green's function in frequency domain 
has no real part. 



5 Schwarzschild-Anti-de Sitter Background 

Since the 3D black hole spacetime is asymptotically anti-de Sitter, one might expect that the late time falloff 
behaviour of any scalar wave in Schwarzschild-anti-de Sitter (SAdS) background has similar behaviour. In 
fact, the situation is quite different from the 3D case, in part because of the different dimensionality and 
in part because it is not possible to solve the wave equation ( pUj ) exactly in this background. Although the 
potential function Ve{r) in this case reads 



VJr) 



^(1-60 + ^ + ^ 



;iio) 
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Figure 16: Decaying exponential behaviour of near the black hole event horizon 



Spinning 3DBH data: Ro = 2R+, RO = 4R+ 
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gure 17: Exponential decay of conformal scalar waves in spinning 3D black hole geometry 
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which is quite similar to (|53|) in the 3D case, we cannot write r as a function of a; in a closed form. That is 
to say given the lapse function 



N{r) 



i_2M ^ \A[^ _Rai±jm 



3 r 3 3 r 

where A < and the black hole radius Rb satisfies the equation N{Rb) = 0, the tortoise coordinate 



'nil 



Rh 



X 



2il + \A\Rb') 



In 



(r-Rb 



r^ + Rbr + Rb' + 3/\A\ 



+ 



V^{2 + \A\Rb 



\A\il + \A\Rb')jA+\A\Rb^ 



arctan 



Ir + Rb 



\V3j4+\A\Rb^ 



TC 

2 



;ii2) 



has no closed form inverse, and so we are unable to write down an explicit expression for the potential V{x). 

For the remainder of this section we set = 1/6 (i.e. we consider a conformal scalar field propagating 
on an SAdS background). 

The most important difference between the 3D black hole background and SAdS background is the shape 
of the potential function V{x). In both cases the potential functions are decreasing in an exponential manner 
toward the horizon. In the 3D case (either spinning or static) this function attains a maximum at a distance 
r = oo [x = 0). When the background is SAdS, this is no longer true because spatial infinity (which is 
still given hj x = 0) is not the place at which V{x) has an absolute maximum (for x G 3?~). As with the 
Schwarzschild black hole (figure the potential function V{x) attains a maximum not far away from the 
event horizon. In the SAdS case, the shape of the potential function V{x) is given in figure 0. Unlike the 
Schwarzschild case, the tortoise coordinate x for the SAdS background is bounded above. Eventually all 
the outgoing waves that leave the black hole region will return towards it due to the boundary condition at 
X = 0. The returning wave will then reflect off of the potential barrier back toward spatial infinity. This 
is completely different behaviour from the 3D case, in which the incoming wave from spatial infinity x = 
continues its journey to the black hole unhindered. 

We can see from figure |T^ that when |A| is small (barriers (a)-(c)), the barrier maximum moves to the left, 
lengthening the travelling time from this maximum to spatial infinity. When / vanishes (case (c) and (C)), 
V{0) = because V^(0) = / (/ + 1) |A|/3 in general. Therefore barriers (a) and (b) have ^(0) 7^ although 
this feature is not apparent on the graph due to the small size of |A|. The barrier height is considerably 
higher than the magnitude of V{x) at x = 0, and this feature becomes more pronounced for large /. This 
causes the scalar wave to bounce back and forth in the region outside the barrier. However part of the scalar 
wave can surmount the barrier (thereby going into the black hole) because the barrier height is still finite. 
However it takes a long time for a significant portion of the wave to enter the black hole. 

We solve the wave equation in this SAdS background numerically. The results of the numerical integration 
of equation (|T0|) using ( |11CI| ) and ( |1 1 l\ j are given in the next few graphs. Figures |T9| and ^ show the falloff 
behaviour of a conformal scalar wave with / = initially located at a distance r = 2 Rb. The computation 
for figure W uses the Dirichlet condition at x = but Neumann boundary condition is employed for figure 
20| . Since the cosmological constant |A| was chosen to be relatively small in both cases, namely |A| = 10"^, 
we can see on both graphs that there is clearly an inverse power-decay behaviour. According to the graph, 
this power-decay rate is roughly t'^ which agrees with the one in Schwarzschild case. However this inverse 
power-decay does not last very long after the return of the outgoing wave from spatial infinity. Both diagrams 
show this returning wavefront. 

For small |A| and nonzero /, the falloff behaviour resembles the case of / = 0. Initially there is a ringing 
effect (due to the quasi-normal modes) followed by inverse power-decay behaviour as shown in figures ^ 
and E2|. 
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Schwarzschild-Anti-de Sitter data: Rb = 2 
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Figure 18: Potential functions V{x) for the SAdS background. The six potentials are generated with the 
parameters (a) A = -10-^/ = 2, (b) A = -10■^^ = 1, (c) A = -10■^^ = 0, (A) A = -10'^, I = 2, (B) 
A = - 10-2, Z = 1, (C) A = - 10-2, 1 = 0. 



SAdS data: Lambda = -0.0001 , Rb = 2 (M = 1 .0001 3), I = 




Figure 19: Scalar wave of / = decays away in SAdS background using Dirichlet condition at x 
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SAdS data: Lambda = -0.0001 , Rb = 2 (M = 1 .0001 3), I = 
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delx = 0.04, delt = 0.03636, Ro = 2Rb, RO = 5Rb 




Figure 20: Scalar wave falloff pattern of I — using Neumann condition spatial infinity 



SAdS data: Lambda = -0.0001 , Rb = 2 (M = 1 .00013), I = 1 




Figure 21: Loglog graph of the decay behaviour in SAdS background using I — 1 and Dirichlet condition 
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SAdS data: Lambda = -0.0001 , Rb = 2 (M = 1 .0001 3), I = 1 
delx = 0.04, delt = 0.03636, Ro = 2Rb, RO = 5Rb 




Figure 22: Similar graph to Figure 21 but using Neumann condition 



Since the finite height of the potential maximum allows the scalar wave to surmount the barrier and leave 
the trapped region which is the exterior outside the barrier maximum, we expect that the peak value of the 
second returning wave is smaller than that of the first wave. However it is unclear from simple inspection of 



figures |T9| to 22 whether or not this is the case. Figures 23 to 26 are the numerical results we obtained using 
the potential functions (A) and (C) in figure We compute the large |A| case with both Dirichlet (figures 
|23| and ^5]) and Neumann and ^6]) boundary conditions. For this case there are more returning waves in 
a reasonable amount of CPU time. From these graphs, we can see that the scalar wave does indeed decrease 
but over a much larger time scale. Figures ^ and ^ indicate that the peak-height has an approximate 
exponential falloff for / = 0. For / > 0, the peak-height has a more complicated behaviour illustrated in 
and |2B. 



figures 
falloff. 



Over long time scales the maximum peak-height has a very mild approximate exponential 



6 Conclusion 

The results of this work indicate clearly that the asymptotics of a given spacetime have a considerable 
influence on the late time falloff behaviour of a scalar wave. The early work of Price showing that waves in 
Schwarzschild background decay away according to an inverse power rate was recently extended by Ching 
et al. who found that the inverse power-law can also be modified by a logarithmic term in some black hole 
configurations. We have been concerned in this paper with the falloff behaviour of conformal scalar waves 
in a background geometry which is asymptotically anti-de Sitter. We found that the waves in uncharged 
static and spinning 3D black hole backgrounds die out at an exponential rate. This conclusion is supported 
by both analytical and numerical computation. For conformal scalar waves in a Schwarzschild-anti-de Sitter 
background our numerical analysis demonstrates that the falloff pattern over short time scales obeys (after 
some ringing) a power-law falloff. However over longer time scales the outgoing wave returns from spatial 
infinity. The falloff of the peak amplitude in this case is neither inverse power nor exponential. However there 
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SAdS data: Lambda = -0.01 , Rb = 2 (M = 1 .01333), I = 



delx = 0.04, delt = 0.03636, Ro = 2Rb, RO = 5Rb 
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Figure 23: Semilog graph of the decay behaviour in SAdS background using I — and Dirichlet condition 
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SAdS data: Lambda = -0.01 , Rb = 2 (M = 1 .01333), I = 



delx = 0.04, delt = 0.03636, Ro = 2Rb, RO = 5Rb 
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Figure 24: Conformal scalar wave decay behaviour using / = and Neumann condition 
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Figure 25: Semilog graph of the decay behaviour using non-zero / and Dirichlet condition 



SAdS data: Lambda = -0.01 , Rb = 2 (M = 1 .00013), I = 2 
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Figure 26: This semilog graph is the Neumann- analogue of Figure 
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Table 1: Exponential falloff rate a which is found by graphical method is always less than the surface gravity 
K in rotating 3D black holes. 



is a very weak exponential decrease in the maximal peak amplitude. This decay pattern is so complicated 
that further investigation will be needed to more precisely determine its dependence on t. 



Previous investigations of mass inflation in a charged spinless 3D black hole |113[ and an uncharged 



spinning 3D black hole |T6[ assumed a power-law falloff rate for the scalar wave. Under such assumptions 
these investigations demonstrated that mass inflation occurred in the corresponding geometries. We have 
shown here that the assumption of power-law falloff in the uncharged static and spinning 3D black hole 
background is not valid and must be replaced by an exponential falloff. 

For a rotating 3D black hole we have shown that a 0-independent scalar field also has late time exponential 
decay. In this particular setup, one can recalculate the mass inflation rate using the correct late time falloff 
rate. It is not difficult to show that the new inflation rate reads m{v) oc exp((K — a) f ), where k is the surface 
gravity at the Cauchy horizon of the black hole. The parameter a comes from the exponential falloff rate 
for the scalar wave (which reads exp(-at)) and is determined (say in the J = case) by setting Iy{ZQ) = 
for the relevant Zq and finding (numerically) the value of u with the maximum negative real part. For 
example when A = -0.1, M = 0.01 and J = 0.001, the surface gravity at the Cauchy horizon is 1.999 but 
the parameter a is only 0.04 according to figure [1^. Table |1| shows some other values of a and k in the 
rotating 3D black hole backgrounds. All these cases are consistent to the claim that the surface gravity k 
is always greater than the exponential falloff rate a. Therefore the attenuation from the falloff effect does 



not halt mass inflation and the basic conclusions in [|16| remain unchanged. Rather the mass inflation rate 



is slowed down relative to spacetimes with power-law falloff. 

We have numerically determined that a is an increasing function of /. Recall that for the approximate 
potential V{x) in (B^), the dominant decay rate of the Green's function is determined by the root ujq of the 



equation Li^/x(ZQ) = in the lower half cu-plane that is closest to the real axis. Since Zq = ^1 + AP/M, 
this root clearly depends upon I. We have computed this dependence and summarized the results in table [| 
, where z/q = -iUo/X- As / increases (i.e. Zq increases), the real part of -Uq also increases, yielding a more 
rapid falloff rate. 

This leads to a qualitative distinction between the mass inflation mechanism in the (static) (2 + 1) and 
(3 + l)-dimensional spacetimes. In a Schwarzschild background, the falloff rate is 1/t^'"*"^ which implies 
larger the moment /, faster the falloff. This yields a mass inflation rate ~ exp(/tf)/f^'+3 [0; the exponential 
growth always surpasses the falloff effect for arbitrarily large /. For a 3D black hole with an inner horizon 
our results imply a falloff rate ~ exp((K — a{l))v). Therefore it is not difficult to imagine a situation in 
which / is so large that the magnitude of the imaginary part of Uq yields an a which exceeds the surface 

*It is interesting to note that the roots with the * in table § are far from precise. With only ten significant figures for i/q, 
the values of I,^^ {Zq) for Zq = 20 and Zq = 30 are of order O(IO^) and O(IO^) respectively. This is the way the modified Bessel 
function behaves for large Zq. In the other cases, the corresponding values of I,y^{Zo) are of order at most 0(10"^). Indeed 
one will need more than 25 significant figures in order to bring 7,^^(30) to the order O(10~^). 
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Table 2: The roots for ly^^Zo) = that are responsible for the dominant exponential falloff 
gravity k which is independent of / at the Cauchy horizon, cutting off mass inflation for those modes. 
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A Solutions in asymptotically flat {D + l)-dimensional back- 
ground 

Given the lapse function (|23|) in (Z}+l)-dimensional spacetime, the time-independent solution of C[ilj{r) ] = 
is given by equation The coefficients aj{r) and Cj(r) can be generated by the next two equations: 



aj+i(r) 



m ( ln|r| )^ CLjir) 

m {D - 2 - a) {D - 1 - a) + I {I + D - 2 - a)] 

27 + 1 

m {D - 2 - a) {D - 1 - a) + {I + a) {I + D - 2)] 



+ 



27 + 1 

m/3 [^{2D -3-2a) + l] 



W(^l + {j + l)a,{ln\r\f aj{r 
1^(27 + 2 + (j + l)a,(ln| 



r\)' aj[r) 



27 + 1 



W[l + ij + l)a,{\n\r\f'^ a^ir) 

mP \n2D-3-2a)-l- D + 2]^^^^ , , ,,x/3-i 

+ ^li^^ ^— ^1^(27 + 2+ j + 1 a, Inr f 

27 + 1 ^ 

"-/^ 1 + (J + 1) «, ( ln|r| a,{r) ) 



aj[r 



27 

m ( ln|r| Y Cj{r) 

m (D - 2 - a) {D - 1 - a) + {I + a) {I + D - 2)] 



(113) 



+ 



+ 



,nlHD-2-a)(D-^l-J) + l(l + D-2-o:)] ^ ^ ^^^^^^ 

ml3[i{2D-3-2a)-l-D + 2 
27 + 1 

m(3 [i{2D-3-2a) + l] 



W^(l + (j + l)a,(ln|r|f ^(r) 

r\f Cjir) 



Vr(l + (j + l)a,(ln|j 



x/3-1 



27+1 



Vr(-27 + (j + l)a,(ln| 



r\f ^ Cj{r-) 
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27+ 

m/5(/5-l)e 



^H/(-27 + (j + l)a,(ln|r|f-'c,(r)) . (114) 



27 + 

The function 1^ used above is defined as 



where n is a real number and /(r) is an integrable function. We evaluate the integral in W in such a way 

that the integration constant is always set to zero. This function has a property that if /(r) is a constant 
and n ^ 1, then W{n, f{r) ) gives a constant. When n is unity but /(r) is still a constant, W{n, f{r) ) is 
proportional to ln|r/7^|, where TZ is some constant of dimension of length. In the case when /(r) is no longer 
constant but a polynomial of ln|r/7^|, the function W gives another polynomial of ln|r/7^|. The index j in 
the two equations runs from zero to infinity and the starting coefficients ao(r) and co(r) are two arbitrary 
constants. It is obvious that when the background is flat, i.e. m = 0, all the %(r) and Cj{r) vanish except 
ao and cq. If m ^ 0, these coefficients are polynomials in ln|r/7^|. 

For the equation £[Ci+i(r) ] = 2dCi{r)/dr with >C[Co(r) ] = 0, the solution of Ci{r) is given by 

air) = r--^-^E^- 
3=0 ' 

Each coefficient c^ (r) can be calculated by the next two equations. 

cHr) - d - 2^27 + 1)! (7 + 0! 
'"^"^^ " '° " 7!(O!(27+l + 0! °' 
c}+\(r) = m {hi\r\f cf^{r) 
2 

+ 



27 + 1 



^ty(-i + (j + l)a,4+i(r)) + (7 + l)W^(-27-l-i + (i + l)«,4+i(r)) 
+ ^[e(^-2-a)(/^-^l^-_.) + (/ + .)(/ + D-2)] , ^ ^ 1) ) 

_..[e(i^-2-.)(Z;-]^.) + U^ + ^-2-c.)] , 

27 + 1 ^ ^ 

+ (, + 1) «, (inM c-(.) ) 

- (- 2 7 - 1 - ^ + (. + 1) «, ( InM )^- c5-(r) ) 

+ ^^^^T^r^ z + (J + 1) a, (ln|r| )^- c5+^(r) ) 
m/3(/3-l)e 



27 + 



^W^(-27-l-i + (j + l)a,(ln|r|)'^-' 4+^(r)) . 
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